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ABSTRACT 

Exoplanetary transmission spectroscopy in the near- infrared using Hubble/NICMOS is currently am- 
biguous because different observational groups claim different results from the same data, depending on 
their analysis methodologies. Spatial scanning with Hubble/ WFC3 provides an opportunity to resolve 
this ambiguity. We here report WFC3 spectroscopy of the giant planets HD 209458b and XO-lb in 
transit, using spatial scanning mode for maximum photon-collecting efficiency. We introduce an analysis 
technique that derives the exoplanetary transmission spectrum without the necessity of explicitly decor- 
relating instrumental effects, and achieves nearly photon-limited precision even at the high flux levels 
collected in spatial scan mode. We validate our method by injecting a synthetic signal into an early 
stage of the analysis, and we recovered it at the correct amplitude. We evaluate our precision using two 
independent error analyses to measure possible red noise. Both methods indicate that our errors are 
within 6% (XO-I) and 26% (HD 209458b) of the photon-limit at a resolving power of X/SX ~ 70. These 
high flux level data allow us to attain a precision better than 0.01% per spectral channel. For XO-lb we 
achieve an average error in (R p /R s ) 2 of 96 ppm, and for HD 209458b our average error is 36 ppm per 
spectral channel over our analyzed range from 1.11-1.65 /im. 

We detect water absorption in both planets, but not other molecules. The absorptions we measure 
are much more subtle than claimed by some previous observers. Both planets exhibit absorption of 
approximately 200 ppm at the water peak near 1.38 fim. Our result for XO-lb contradicts the much 
larger absorption derived from NICMOS spectroscopy. The weak water absorption we measure for 
HD 209458b is reminiscent of the weakness of sodium absorption in the first transmission spectroscopy of 
an exoplanet atmosphere by Charbonneau et al. (2002). Model atmospheres having uniformly-distributed 
extra opacity of 0.012 cm 2 g _1 account approximately for both our water measurement and the sodium 
absorption in this planet. Our results for HD 209458b support the picture advocated by Pont et al. (2012) 
in which weak molecular absorptions are superposed on a transmission spectrum that is dominated by 
continuous opacity due to haze and/or dust. However, the extra opacity needed for HD 209458b is grayer 
than for HD 189733b, with a weaker Rayleigh component. 

Subject headings: stars: planetary systems - transits - techniques: photometric - techniques: 
spectroscopic 
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1. INTRODUCTION 

P rom the first discovery of transiting extrasol ar plan- 
ets (jCharbonneau et all l2000t iHenrv etaH I2000H trans- 
mission spectroscopy was anticipated a s a potential tech- 
nique to probe their atmospheres (Sea ger fc Sasselcvl 
2000). Indeed, transmission spectroscopy was used to 
make the first detection o f an exoplanetary atmosphere 
(|Charbonneau et all [2002), via optical sodium absorp- 
tion observed using the Hubble Space Telescope (HST), 
and sodium and potassium absorpti on measurements are 
now possible from the ground (e.g., iRedfield et all [2008; 



201J) - Using HST 



ISnellen et all 120091: ISing et all 120111 ._ 
data, Barman! 1 20071) identified water abs orption near 
1 /xm i n the giant exoplanet HD 209458b, and lDesert et all 
(2008) searched for evidence of TiO/VO absorption in 
that planet. Expanding HST trans mission spectroscopy 
to longer infrared (IR) wavelengths, I Swain et all (|2P08a) 
obtained results indicating water and methane absorption 
near 2 /xm in the gi ant exoplanet HD 189733b. Similarly, 
iTinetti et all (|2010f ) derived water and carbon dioxide ab- 
sorption near 1.4 /xm during the transit of XO-lb. 

Successful transmission spectroscopy of giant exoplan- 
etary atmospheres is a crucial first step toward even- 
tual spectroscopy of a nearby habitable sup er-Earth using 
the J ames Webb Space Telescope (JWST, iDeming et all 
2009). However, the IR transmissio n spectroscopy using 
HST/NICMOS (ISwain et all l2008afc ITinetti et all [2010ft 
has been challenged bv lGibson et al.l (|201lh who empha- 
size that the reported absorption features are sensitive to 
corrections for instr umental systematic errors. Moreover, 
iGibson et all (l201lh argue that corrections for instrumen- 
tal error cannot be made using simple linear decorrela- 
tions because the 'instrument m o del' i s not sufficiently 
well understood. ICrouzet et al.l (|2012f ) concluded that 
NICMOS instrumental signatures remain comparable with 
the expected amplitude of molecular signatures, even af- 
ter a decorrelation analysis. Since the HST/NICMOS 
instrument is no longer in operation, repeating the ob- 
servations to clarify this situation is impossible. Fortu- 
nately, transmission spectroscopy from 1.1 to 1.7 /xm - 
largely overlapping the NICMOS G141 grism region at 1.2- 
1.8 /xm - i s possible using Wide Field Camera 3 (WFC3) 
on HST (jBerta et al.l 120121 IGibson et all [2012I ). More- 
over, the WFC3 detector is known to exhibit a more uni- 
form intra-pixel se nsitivity response than does NICMOS 
(|McCulloughll2008l) . giving reason to expect that WFC3 
observations may be less affected by instrumental signa- 
tures. Also, WFC3 can now be operat ed in a spatial scan 
mode (|McCullough Sz MacKentvll2012l ) wherein the target 
star is trailed during each exposure by telescope motion 
perpendicular to the direction of dispersion. Exoplanet 
host stars are often bright, and the spatial scan allows the 
longest practical exposures for bright stars without satu- 
rating the detector, greatly increasing the overall efficiency 
of the observations. 

In this paper we report WFC3 transmission spec- 
troscopy for two exoplanets observed in our Large (115 or- 
bit) HST Cycle-18 program. By coincidence, these planets 
(XO-lb and HD 209458b) were both scheduled for obser- 
vation late in our program, permitting us to acquire the 
spectra in the newly-developed spatial scan mode. One of 
these planets (XO-lb) is the same as observed in trans- 



mission using NIC MOS spectroscopy (ITinetti et all l2010l: 
IGibson et alj|201ll ICrouzet et all |2012|). In addition to 
the great photon-collecting efficiency provided by spatial 
scan mode, we have achieved some new insights in the 
analysis of W FC3 data, beyond t he valuable methodology 
introduced bv lBerta et ail (|2012l ). 

We here report robust exoplanetary transmission spec- 
tra using WFC3 in spatial scan mode. We achieve a level 
of precision closely approaching the limit imposed by pho- 
ton statistics, even for these large exposure levels collected 
in spatial scan mode. Moreover, our analysis methodol- 
ogy requires no explicit decorrelation using an 'instrument 
model', nor does it require that the pattern of systematic 
errors be consistent between orbits, or that we omit the 
first orbit per visit from our analysis. 

Sec. 2 describes the circumstances of our observations, 
and Sec. 3 discusses the initial calibration of the data, 
including a brief discussion concerning the nature of in- 
strumental signatures produced by WFC3. In Sec. 4 we 
describe new methodology that we have used to extract 
the transmission spectra of the planets, and Sec. 5 gives 
our results and relates them to previous work. Sec. 6 inter- 
prets our results using model atmospheres for the planets, 
and Sec. 7 summarizes and comments on future possibili- 
ties. 

2. OBSERVATIONS 

Observations of HD 209458b and XO-lb used WFC3 
with the G141 grism, providing wavelength coverage from 
1.0 - to 1.7 /xm in first order. Each star was observed dur- 
ing a single visit comprising 5 consecutive orbits. The 
observational sequence in each visit began with an undis- 
persed image of the star using the F139M filter (central 
wavelength of 1390 nm, and FWHM = 32 nm). The filter 
choice for this image is not crucial, because its purpose is 
merely to establish the position of the undispersed stel- 
lar image, used in wavelength calibration. Following the 
undispersed exposure, we obtained a sequence of exposures 
using the grism, with the star scanned perpendicular to the 
dispersion, under control of the HST fine-guidance system. 

The grism exposures used subarray readouts of the de- 
tector to maximize efficiency. Information on the number 
of exposures, the duration of each sequence, subarray size, 
and range of planetary orbital phases are given in Table 1. 
Table 1 also includes the spatial scan parameters and av- 
erage exposure level in the spectral images, because this 
information may be useful to subsequent observers. 

3. INITIAL DATA PROCESSING 

WFC3 grism spectroscopy samples the detector 'up the 
ramp', i.e., reading each pixel non-destructively multiple 
times. We process these data by a method described in 
the Appendix. After the initial processing, the spectral 
frames have wavelength in one dimension, with the spatial 
scan in the orthogonal direction. 

Examples of the 2-D spectral frames are shown in Fig- 
ure 1. One characteristic of spatial scan mode is that the 
rate of telescope motion is not perfectly uniform, but varies 
slightly with time due to jitter in the control by the fine 
guidance system. Evidence of this variation is seen on 
the rightmost image of Figure 1, that shows the difference 
between two consecutive scans of HD 209458b. The vari- 



Exoplanetary Transmission Spectroscopy 



3 



able scan rate results in variable flux as a function of row 
number, typically varying by ±1% as shown in the plot on 
Figure 1. Fortunately, because the scan is perpendicular to 
the direction of dispersion, this phenomenon does not con- 
tribute significant noise to our analysis, but it does affect 
how we process the data to discriminate against hot pixel 
and energetic particle hits, as the Appendix describes. 

3.f . Wavelength Calibration and Flat-fielding 

Following the initial data processing described above, 
we apply wavelength and flat-field calibrations. Each 
visit includes an undispersed image of the star, and the 
wavelengths in the grism spectrum are reliably fixed rel- 
ative to the position of the undispersed image. We cal- 
culate thos e wavelengths using coe fficients recommended 
by ST ScI dKuntschner et al] |2009D. with some modifica- 
tions. iWilkins et al.l (12013ft found that the original rec- 
ommended coefficients did not produce optimal agree- 
ment between the overall profile of the grism response 
in observed spectra when compared to the known profile 
of that response, and also did not give exact agreement 
with the known wa velengths of stellar absorption lines. 
IWilkins et al.l ()2013l ) varied the coefficients by up to 10% 
by trial- and-error, to achieve o ptimal agreement. We have 
used the Wil kins et al.l (|2013h coefficients, and we obtain 
good consistency of the grism response, and good account- 
ing of stellar absorption lines (e.g., Paschen-beta, see be- 
low). 

In addition to wavelength calibration, we apply the 
wavelength-dependent flat field calibration as recom- 
mended by STScI. Note that flat-fielding is not included in 
STScI pipeline processing, even for the files of non-spatial- 
scan data, because it is a function of where each target's 
spectrum happens to fall on the detector, and must be 
done at the user stage of data analysis. To the extent 
that the stellar spectra were fixed on the detector during 
each visit, and jitter and drift in wavelength were negligi- 
ble, flat-fielding would not be needed in our analysis. In 
that (ideal) case, we would not apply the flat-fielding step, 
since it has been our experience that the more the data are 
processed, the more difficult it becomes to achieve photon- 
limited results. However, wavelength jitter in spatial scan 
mode can be larger than in non-scanned observations, so 
the flat-fielding step of the analysis is prudent. However, 
we have repeated the entire analysis of this paper without 
the flat-fielding step, and we find consistent results in the 
two cases. 

3.2. Instrumental Signatures 

Our analysis is designed to be insensitive to effects 
caused by the instrument and detector. To understand 
how we minimize such effects, we must briefly discuss 
what is known about instrumental signatures in WFC3 
GI4I grism data. Figure 2 shows a portion of the transit 
of WASP-I8b, which was observed in our program using 
non-scanned mode, and which clearly illustrates this ef- 
fect (our spatial scan data show it less clearly). The dom- 
inant systematic error is an increase in intensity during 
each group of exposures that are obtained between buffer 
dumps. This pattern is shaped like a T', and is slightly 
reminiscent of a fish hook. It may be physically similar to 
the ramp effect seen prominently in Spitzer 8 /im data, but 



the WFC3 time scale is shorter and we cannot be certain 
that it has the same physical cause as the Spitzer ramp. 
We therefore use different terminology, and refer to this 
WFC3 phenomenon as the 'hook'. 

The amplitude of the hook has been studie d as a func- 
tion of exposure level bv IWilkins et al.l (12013D . They find 
that the hook is, on average, zero when the exposure level 
per frame is less than about 30,000 electrons per pixel. 
Above that exposure level, the amplitude of the hook 
increases with greater exposure levels, albeit with rela- 
tively large scatter in the relation. Our exposure levels 
are about 40,000 electro ns per pixel (Table 1). Based on 
the IWilkins et al.l (|2013h results, we expect the hook to be 
weakly detectable in our data. 

4. DERIVATION OF THE EXOPLANETARY SPECTRA 

The principle of our analysis is that we first define and 
fit the transit as observed in the integral intensity of the 
grism spectrum over wavelength. Our subsequent analysis 
removes this 'white light' transit, and solves for differential 
transit depths as a function of wavelength. We construct 
the exoplanetary transmission spectrum by adding these 
differential amplitudes to the white light amplitude to pro- 
duce the transit depth versus w avelength. There i s prece - 
dent for this methodology, from iRichardson et al.l (|2007j ). 
In what follows, we describe the steps of our analysis in 
detail, but we also provide a concise summary in Sec. 4.7. 

Our work has benefitted from experi ence and analy- 
ses of non-scanned d a ta in our prog r am dLine et al.ll2013l: 
i Mandell et ail 120131: iRanian et all 120131: IWilkins et al l 
| 2013h. as wel l as th e analyses bv lBerta et al.l ( 20121 ) and 
Gib son et al.l (|2012f) . Some of the lessons learned in those 
analyses are: that the hook is usually common-mode to 
different wavelengths and will cancel in an appropriate ra- 
tio. Also, derivation of the exoplanetary spectrum at the 
edges of the bandpass can be problematic for two reasons. 
First, the grism spectrum exhibits jitter in wavelength that 
especially affects the strongly-sloped edges of the band and 
must be corrected as part of the analysis. Second, the in- 
tensities decrease at the edges of the grism response, so the 
amplitude of the hook for pixels sensing those wavelengths 
will be less, and potentially not common-mode with the 
remainder of the spectrum. 

Following the application of wavelength and flat-field 
calibration to the 2-D spectral frames, we sum each spec- 
tral frame spatially (i.e., along columns) to derive 1-D 
spectra. Mindful of the experience related above, we per- 
form the spatial sum over a range of rows slightly less than 
the full spatial extent of the scan. In this manner, we avoid 
including lower intensities that occur at the edges of the 
scan. We sum in wavelength over a range that exceeds 
the region of significant grism response, and we restrict 
the wavelength range (and thus intensity range) in a later 
step of the analysis (see below). We include the first or- 
bit of each visit in our analysis, but as a precaution we 
omit the first 5 frames of the first orbits, where the hook 
is strongest. 

4.1. The White Light Transit 

Figure 3 shows the white light photometry for our ob- 
served transits of each planet. We have not attempted to 
correct for the hook, using a divide-out-of-transit (divide- 
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oot) procedure as described by iBerta et al.l ()2012l ). The 
hook in these spatial scan data is weak and inconsistent 
(Figure 3), and good fits are possible without making cor- 
rections. To fit the white light transit of HD 209458, we 
ge nerate theoretical tran sit curves using the formulation 
of lMandel fc Agoll t ()02). using non-linear limb darkening 



coefficients ( Claretl feOOOf) . calculated from the emergent 



intensities of a 3-D hydrodyn amic model of the stellar at- 
mosphere (|Havek et al.l l2012h . We integrated the model 
intensities over our specific WFC3 bandpass, and fitted 
the results to the nonlinear law to obtain the coefficients. 
To quantify the degree to which the results for HD 209458 
depend on the treatment of limb darke ning, we also do the 
fits u sing a linear limb-darkening law (jClaret fc Bloemenl 
[20ll . We find good agreement in results between the lin- 
ear and nonlinear limb darkening treatment for HD 209458 
(see below). Nonlinear coefficients for XO-1 would require 
more interpolation in the 3-D model atmosphere grids than 
we prefer, and our observed noise level is higher for XO-1 
than for HD 209458. We therefore use linear limb darken- 
ing for XO-1. For both HD 209458 and XO-1, we interpo- 
late the linear coefficients based on values from adjacent 
bandpasses (J and H). 

For the linear limb darkening coefficients, a range of val- 
ues are available that vary (typically by ±0.03) depending 
on whether they are calculated from PHOENIX versus AT- 
LAS models, and depending on minor parameters such as 
microturbulence, and on the method used to fit the model 
atmosphere to the linear law. Our adopted values for lin- 
ear coeffic ients correspond to t he mi ddle of the range tab- 
ulated by IClaret fc Bloemenl (|201l[ ) . A similar range of 
non-linear coefficients is not available, so we use the varia- 
tion of the linear coefficient to estimate the impact of limb 
darkening uncertainties. 

For fitting the HD 209458b transit, we adopt the o rbital 
parameters (P, i, a/R s ) from lKnutson et all ()2007al ). with 
linear limb darkening coefficient 0.28. For the non-linear 
case we use coefficients calculated for our specific band- 
pass, as noted above. Our WFC3 observations show a 
visit-long downward trend. Omitting the first orbit, we 
are able to fit the remaining data using a linear baseline 
plus the transit. (Note, however, that we do not omit the 
first orbit from the spectrum derivation discu ssed below.) 
For X O-lb we use the orbital parameters from lBurke et al.l 
(20ld[), with linear limb darkening coefficient of 0.31. We 
include the first orbit for XO-1. For each planet, with 
the limb darkening and orbital parameters fixed, we vary 
Rp/R s and the time of transit center (to account for im- 
precision in the ephemeris). The resulting best- fit transit 
curves are overplotted (blue lines) on Figure 3. 

To estimate the errors on our derived val ues, we use a 
resid ual permutation ('prayer-bead') method ()Gillon et al.l 
2007). Also, we manually vary the lin ear limb darkening 
coeffic ient over the range tabulated bv lClaret fc Bloemenl 
(|201lh , and note the impact on R p /R s . Our adopted errors 
are the quadrature sum of the variation as a function of lin- 
ear limb-darkening coefficient, and the prayer-bead errors. 
For HD 209458b we derive R p /R s = 0.1209 ± 0.0005 using 
non-linear limb darkening, and R p /R s — 0.1214 ± 0.0005 
using linear limb darkening. Within th e errors, our re- 
sults a gree with other IR transit results. ICrossfield et al.l 
(|2012[) derived R p /R s = 0.1218 ± 0.0014 at 24 /im from 



several Spitzer transits. For XO-lb, we derive R p /R s = 
0.1328 ± 0.0006, that agrees (within the errors) with a 
transit observed using NICMOS by iBurke et all (|2010h 
(Rp/Rs = 0.1320 ± 0.0005). Our retrieved transit times 
are close (tens of seconds) to the predictions using the or- 
bital parameters cited above. Our results for radius ratio 
and transit time are summarized in Table 2. 

4.2. Wavelength Shifts 

Upon summing the 2-D frames in the scan direction to 
produce 1-D spectra, we find that these spectra are not 
coincident in wavelength as judged by the displacement 
of the grism response curve. Variations up to ±l-cohmm 
(i.e., in the wavelength direction) occur over the span of 
each visit. These variations are almost two orders of mag- 
nitude larger than similar variations (±0.02-columns) seen 
in non-scanned data. That is reasonable, since it should 
be much easier for the fine guidance system to hold a fixed 
position in wavelength when it does not have to scan spa- 
tially. 

A high quality analysis of spatially scanned data must 
account for the wavelength shifts. We proceed as follows. 
First, we construct a template spectrum by averaging spec- 
tra that occur before first contact or after fourth contact, 
by up to one hour. By restricting the temporal range of 
these out-of-transit spectra to within one hour of the con- 
tacts, we minimize potential broadening of the template 
due to the wavelength shifts themselves. After forming 
the template spectrum we shift it in wavelength, and fit it 
to all of the individual grism spectra in the visit. For each 
grism spectrum, we step through a large range of wave- 
length shifts of the template, using 0.001-pixel increments, 
to assure that the best-fit shift is identified. At each shift 
value, the best trial fit is achieved by a linear stretch of the 
template spectrum in intensity (via linear least-squares). 
After trying the full range of possible wavelength shifts, 
we pick the best wavelength shift and linear stretch factor 
based on the minimum \ 2 ■ Shifting the template spectrum 
requires resampling it by interpolation. We use the IDL 
routine INTERPOLATE, with the cubic keyword set to 
-0.5, this being the best approximation to ideal interpola- 
tion using a sine function. Because we shift the template 
to match individual spectra, re-sampling of each individ- 
ual spectrum is avoided, further minimizing the potential 
for adding noise via the re-sampling process. 

Upon fitting the template spectrum to a given individ- 
ual spectrum, we subtract them and form residuals R t \, 
where t indexes time (i.e. what individual spectrum), and 
A indexes wavelength. This subtraction removes the small 
amount of sky background that survives the process de- 
scribed in the Appendix. Also, we normalize the Rt\ by 
dividing by the template spectrum. An illustration of the 
match between an individual spectrum and the template 
is shown in Figure 4, including the residuals in the lower 
panel. At each A, the time series Rt\ contain the differen- 
tial transits that we seek. To facilitate cancellation of the 
hook and potentially other systematic errors, we restrict 
the intensity range of our analysis - hence the wavelength 
range - to wavelengths whose intensity lies above the half- 
power points of the grism response (see dashed lines on 
Figure 4). 

The normalization of the R t \ is key to our analysis, be- 
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cause it helps to cancel common- mode systematic errors. 
That cancellation is conceptually equivalent to dividing 
the intensity in the grism spectrum at a given wavelength 
by the integral of each grism spectrum over wavelength. 
Thus, o ur methodology is similar to the divide-oot method 
used bv lBerta et al.1 (|2012t) , but is (arguably) more general. 
The divide-oot method relies on the pattern of systematic 
error being consistent in time, whereas we here require 
that it be common-mode in wavelength. 

To see why our procedure described above is equivalent 
to a wavelength ratio, consider the following. We use a 
single template spectrum per planet for each visit. We 
stretch that template spectrum in the fitting process, but 
the stretch is a linear factor common to all wavelengths. 
Therefore, the ratio of the template at a given wavelength 
to its integral over wavelength is constant. Moreover, the 
fitting process guarantees that the wavelength integral of 
the stretched template will be closely equal to the wave- 
length integral of the individual spectrum being fit. Hence, 
apart from a constant factor, dividing an individual resid- 
ual by the value of the template spectrum at that wave- 
length is conceptually the same as dividing the individual 
spectrum at that wavelength by its integral over wave- 
length. 

Although we have described our method as being equiv- 
alent to a ratio-of- wavelengths, the illumination level of the 
various pixels is probably a more relevant physical variable 
than wavelength per se. Since our analysis is restricted to 
intensities not greatly below the peak of the grism response 
(Figure 4), the intensities in the data covered by our anal- 
ysis tend to be restricted to a limited range, and this is 
probably the dominant factor in cancellation of systematic 
errors. 

Note that our analysis procedure as described above 
(and further below) does not involve any explicit decor- 
relation versus instrument or telescope parameters (e.g., 
tilt of the spectrum, Hubble orbital phase, d etector tem- 
perat ure, etc.). Like the divide-oot method (|Berta et al.l 
I2012T ). we rely on cancellation of common- mode system- 
atic errors by operating only on the data themselves, using 
simple linear procedures. 

4.3. Undersampling 

Initial correction of wavelength shifts using the above 
procedure showed discordant results at some wavelengths, 
characterized by strong slopes and even non-linear tem- 
poral trends in the R t \ values as a function of t. The 
most discordant results occurred near strong stellar lines 
such as Paschen-beta (1.28 /xm). We initially suspected 
interpolation errors in the shift-and-fit process, but care- 
ful inspection of the uninterpolated spectra revealed that 
the shapes of the stellar lines were changing as a func- 
tion of the wavelength shift, due to undersampling of the 
grism resolution by the pixel grid. The FWHM resolution 
of the G141 grism at 1.28 /xm equals 2.3 pixels. This is 
insufficient to eliminate changes in the pixel-sampled line 
shape as a function of wavelength shift. Figure 5 shows 
the Paschen-beta line in two spectra separated by about 3 
hours in the visit for HD 209458. The change in line shape 
is obvious. This line shape change is not mere noise, since 
it is consistent over many spectra, and the pixel-sampled 
line shape changes gradually and smoothly as a function 



of wavelength shift. 

WFC3 sampling of 2.3 pixels per spectral resolution 
nominally complies with th e criterion of t he Nyquist- 
Shannon sampling theorem (|Shannonl Il949f ). However, 
the grism spectral response can violate the premise of the 
theorem in the sense that its Fourier decomposition may 
contain components at spatial frequencies higher than the 
nominal resolution. So the undersampling we infer here is 
not surprising. 

Changes in the shape of pixel-sampled stellar lines dur- 
ing a transit will cause noise that cannot be removed using 
any simple divide-oot or ratio procedure. Moreover, we 
cannot change the dimensions of the pixel spacing. Our 
solution is to force adequate sampling of the spectral res- 
olution, by degrading the resolution post-detection. Prior 
to the analysis described above (i.e., before forming and 
using the template spectrum), we convolve each 1-D spec- 
trum with a Gaussian kernel having FWHM = 4 pixels. 
We varied the width of the kernel to determine the best 
compromise between supression of undersampling errors 
and degradation of the spectral resolution. We apply the 
convolution to the template spectrum as well as to in- 
dividual spectra. Figure 4 has been convolved with our 
adopted kernel. Because of the linearity of our analysis, 
we arguably could achieve similar results by fitting the 
template spectrum in the presence of undersampling er- 
rors and averaging the resultant exoplanetary spectrum 
over wavelength in a subsequent step of the analysis. We 
elect to smooth the grism spectra at an early stage of the 
analysis because that gives us more insight into the nature 
of the errors when deriving differential transit amplitudes. 

4.4. Differential Transits 

The wavelength dependence of the transit depths is con- 
tained in the Rt\ residuals. Note that the smoothing pro- 
cedure described above introduces autocorrelation in the 
residuals as a function of A (apparent on Figure 4), but 
not as a function of t. At each A, we fit a scaled transit 
curve to the Rt\, with a linear baseline (i.e., we allow for 
a linear trend with time). The shape of the transit curve 
is constrained to be the same as the white light (Figure 3) 
curve for each planet, except that we include a correction 
for the wavelength dependence of limb darkening. We al- 
low the amplitude of the fitted curve to vary, since that is 
essential to deriving the exoplanetary spectrum. The best- 
fit amplitude and baseline slope are found simultaneously 
via linear regression. The slopes retrieved from the regres- 
sion are small, and have little impact on our r esults . Note 
that more sophisticated (e.g., Markov-Chain. iFordl 12005) 
techniques would be superfluous in this situation, since we 
are not concerned with errors introduced by correlations 
between parameters, etc. Nor do we need to consider un- 
certainties on priors like the planet's orbital parameters, 
because the relevant priors are already known to high pre- 
cision and are not dominant in our analysis. Instead, the 
dominant problem is simply to find the best-fit diffcrental 
amplitude and baseline in the presence of noise. 

Although we calculate the differential transit amplitudes 
by fitting a transit curve with a linear baseline, we also 
checked the results using a much simpler procedure. Divid- 
ing the R t \ at each A into an in-transit and out-of-transit 
group, we subtract the average of the out-of-transit resid- 
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uals from the average of the in-transit residuals at each 
A. This simple in-minus-out procedure yields transmission 
spectra that are very similar to the more rigorous method 
of fitting a transit curve (fitting the curve accounts for the 
ingress and egress portions correctly, and it permits us to 
correct for the wavelength dependence of limb darkening). 

There has been considerable discussion in the liter- 
ature concerning method ologies to derive exoplanetary 
transm ission spectra (e.g. . [Gibson et al1l201ll : ISwain et al.l 
2008 jl), including some quite sophistica ted techniques 
( Gib son et al J 12012: IWddmann et al.ll20H . While we re- 
spect the power of sophisticated analyses, we advocate the 
virtue of making the signal visible to the eye using the sim- 
plest linear processes. To that end, we present Figure 6, 
that shows the differential transit data for HD 209458b, 
averaged over two wavelength regions - one region of high 
water opacity (1.36-1.44 /xm) and a second region of low 
water opacity (1.27-1.30 /xm). These transit curves are av- 
eraged over multiple A values (17 and 6 columns of the de- 
tector, respectively), whereas our actual differential transit 
analysis solves for differential amplitudes at individual A 
values. However, the averaging used for Figure 6 is lin- 
ear (fit to the average is the average of the fits), and the 
modest averaging helps to make the signal obvious by eye. 
Note that the region of high water opacity shows an ob- 
viously deeper-than-average transit, and the low opacity 
region exhibits a negative differential amplitude, i.e. a 
less-than-average transit depth. 

4.5. Verification of Sensitivity 

Anticipating our results (Sec. 5), we find exoplanetary 
water absorptions that are of significantly smaller ampli- 
tude than previous investigators claim for the same plan- 
ets. We therefore verified the sensitivity of our analysis by 
numerically injecting a synthetic signal into our data, and 
we recovered it with the correct amplitude. Specifically, we 
added a synthetic transit of amplitude 500 ppm, occurring 
only in 10 columns of the detector spanning wavelengths 
1.225-1.272 /xm, to the HD 209458b data. We added this 
synthetic signal immediately after the stage of producing 
the scanned data frames (Eq. 1 in the Appendix). Our 
analysis retrieved this signal at the full injected ampli- 
tude (not illustrated here), with the expected roll-off at 
the edges of the simulated sharp band due to the smooth- 
ing used in our analysis. We conclude that our analysis 
does not numerically attenuate exoplanetary transmission 
signals to any significant extent. 

4.6. Errors 

We estimate the errors associated with our differential 
transit amplitudes using two techniques. First, we calcu- 
late the standard deviation of single points in each differ- 
ential transit curve, after the best-fit differential transit is 
removed (i.e., in the residuals). We denote this value as 
<ti . Then, we bin the residuals of each transit curve over N 
points, varying N up to half the number of observed points 
in time, and we calculate the standard deviation of each 
binned set, denoted er/v. For Poisson errors in the absence 
of red noise, we expect: 

logiajf) = /o 5 (<7i) - 0.5 log(N) (1) 



The slope of the observed log(ai) versus log(apf) rela- 
tion is uncertain at a single wavelength due to the paucity 
of points at high-N (i.e., few large bins). For better statis- 
tics, we accumulate the a values over all wavelengths, and 
show the dependence of ctn versus log(N) on Figure 7 
for both planets. The blue lines on Figure 7 represent 
the errors expected in the limit of photon statistical noise, 
that decreases proportional to —0.5log(N). These lines 
are calculated a priori from the number of electrons in 
the spectrum, and their close accord with the measured 
scatter implies that the errors of our analysis are close to 
photon- limited. Our cumulative o~(N) values are in good 
agreement with the expected photon noise and the slope 
of -0.5 in Eq.(2). One seeming difference on Figure 7 is 
that both planets exhibit some cr/v points that scatter well 
below the photon noise limit at large AT. This occurs be- 
cause the differential transit fitting acts as a high-pass fil- 
ter. Even if the differential transit amplitude is zero, the 
linear regression will often find a non-zero amplitude due 
to noise at large bin sizes. The regressions will therefore 
tend to remove low frequency noise as a by-product of de- 
riving the differential transit amplitudes. 

Based on Figure 7, we calculate the error associated with 
each differential transit amplitude as being the quadrature 
sum of the in- and out-of-transit levels in the Rt\ values, 
calculating the errors on these levels from Eq.(2) . We also 
used the prayer-bead method (|Gillon et alJ uOO?) to calcu- 
late the error on each differential transit amplitude. This 
also indicated close agreement with the photon limit, but 
the precision of the prayer-bead error estimate at a sin- 
gle wavelength is limited by the relatively small number 
of possible permutations. Therefore we adopted our final 
errors using the following procedure. For each planet, we 
calculated the average ratio of the prayer bead to photon 
errors, and the scatter about this average. We multiply the 
photon errors by this average ratio to obtain our adopted 
errors for most wavelength bins. However, a few wave- 
lengths exhibit prayer bead errors more than 3tr greater 
than the average error level. For these points, we use the 
prayer bead error estimate for that specific wavelength. 

Our derived transmission spectra for HD 209458b and 
XO-lb are shown on Figure 8. Our average error bar for 
the spectrum of HD 209458b is 36 ppm, which is 1.26 times 
the photon noise. The largest error bar for our 28 wave- 
length channels is 51 ppm. For XO-lb, our average errors 
are 96 ppm, which is 1.06 times the photon noise. The 
largest XO-lb error bar is 110 ppm. 

4.7. Summary of Our Spectral Derivation Methodology 
To summarize our method as described above: 

• From 2-D spectral images that are fiat-fielded and 
wavelength-calibrated, we make 1-D grism spec- 
tra. We sum the 2-D spectral images over a range 
slightly less than their height, to utilize pixels hav- 
ing similar exposure levels, to the maximum pos- 
sible degree. We similarly restrict our analysis to 
wavelengths well above the half-intensity points on 
the grism sensitivity function, also to use pixels 
with similar exposure levels as much as possible. 

• We integrate the grism spectra over wavelength 
within our adopted wavelength range, and con- 
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struct a band-integrated transit curve. We fit to 
this transit curve to obtain the white-light tran- 
sit depth (Rp/R^). We save the white-light transit 
depth to use below. 

• We smooth the grism spectra using a Gaussian ker- 
nel with a FWHM = 4 pixels. This reduces the 
effect of undersampling. We construct a template 
spectrum from the out-of-transit smoothed spectra, 
and we shift it in wavelength, and scale it linearly 
in intensity, to match each individual grism spec- 
trum, choosing the best shift and scale factors using 
linear least-squares. We subtract the shifted and 
scaled template to form residuals, and normalize 
the residuals by dividing by the template spectrum. 
This procedure removes the white-light transit, but 
preserves the wavelength variation in transit depth. 

• At each wavelength, we fit a transit curve to the 
residuals as a function of time, accounting for the 
wavelength dependence of stellar limb darkening. 
We add the amplitude of this transit curve (a 'dif- 
ferential amplitude') to the depth of the white-light 
transit from above. We then co-add the results in 
groups of 4 wavelengths (columns on the detector) 
to match the smoothing described above. The re- 
sult is the exoplanetary transmission spectrum. 

• We determine errors using a residual-permutation 
method, comparing those to errors calculated by 
binning the residuals over increasing time intervals 
(to verify an inverse square-root dependence), and 
by comparing to an ab initio estimate of the photon 
noise. 

• We verify the sensitivity of the method to assure 
that it does not numerically attenuate the exoplan- 
etary spectrum. We inject numerically an artifi- 
cial spectrum into the data at the earliest practical 
stage of the analysis, and we recover it at the cor- 
rect amplitude. 

5. RESULTS FOR TRANSMISSION SPECTRA 
5.1. Comparison to Expectations from Spitzer 

We here illustrate our observed spectra, and immedi- 
ately compare them to models that are consistent with 
Spitzer emergent intensity observations of these planets. 
The strong similarity between observed and modeled spec- 
tra reinforces our conclusion (Sec. 6.1) that we are observ- 
ing real exoplanetary absorption. In Sec. 6, we explore 
comparisons with model atmospheres in more depth. 

Figure 8 overplots modeled transmission spectra for each 
planet in comparison to our observed result s. The modeled 
spectra were calculated by Adam Bu rrows (|Burrows et al.l 
120011 120101: iHowe fc Burrows! I2012T ) based on combining 
day- and night-side model atmospheres that are consistent 
with Spi tzer secondary eclipse observations for the se two 
planets ([Burrows et all 120071: [Machalek et alll2008D . The 
day and night-side model atmospheres were combined by 
equalizing their basal pressures to join them at the ter- 
minator of the planet. The transmittance spectrum used 
for Figure 8 represents a line of sight that passes through 



both the day- and night-side models. We fit the mod- 
eled absorptions to the data by scaling them in ampli- 
tude, and offsetting them slightly in overall radius, but 
not changing the modeled shape. The fitted amplitude of 
the HD 209458b absorption is 0.57 of the model predic- 
tion, and for XO-lb the fitted amplitude is 0.84 of the 
model prediction. These factors are physically reasonable, 
as we discuss below, and the correspondence between the 
observed and modeled shape of the absorptions is clear. 

Our analysis uses simple procedures without recourse to 
an instrument model, and the smoothing we implement is 
motivated by an effect that we understand physically (the 
undersampling) . Our errors are close to photon-limited as 
verified by the inverse square-root dependence when bin- 
ing, and by the comparison with the prayer-bead errors. 
We are therefore certain that Figure 8 represents the real 
astrophysical absorption spectrum, especially in the case 
of HD 209458b, where the amplitude of the absorption (200 
ppm) is more than 5 times the average error per point, and 
the absorption is sampled by many observed points. 

In the case of XO-lb, the detection is less visually obvi- 
ous than for HD 209458b, but is still robust. If we fit a flat 
line (the no absorption case) to the XO-lb observations, 
the x 2 is 64.6 for 28 degrees of freedom. That rejects 
the flat line at greater than 99% confidence. Moreover, 
the total context of the observations, including the simi- 
larity to both the model prediction and the HD 209458b 
observations, allows us to conclude that real astrophysical 
absorption is observed in XO-lb as well as in HD 209458b. 

In Sec. 6.1 we consider whether stellar activity could 
contribute significantly to our derived exoplanetary spec- 
tra, or whether true wavelength-dependent absorption in 
the planetary atmospheres is dominant, and we conclude 
the latter. A more elaborated comparison with planetary 
models is presented in Sec. 6. 

5.2. Comparison to Previous Observational Results 

The absorptions we detect are considerably weaker 
than claimed by several previous investigations. The 
clearest discrepancy is for XO-lb, as illustrated on Fig- 
ure 9. This repeats our XO-lb spectrum, ove r laid o n the 
same plot as the re sults from iTinetti et al.l (|2010f ) and 
ICrouzet et all (l2012[). T he large water absorption derived 
by ITinetti et all "! 20101) is inconsistent with our results: 
such a large signal w ould be obviou s even at the Figure 4 
stage of our analysis. ICrouzet et al.l (|2012l ) concluded that 
NICMOS instrumental signatures remain comparable with 
the expected amplitude of molecular features , even after a 
decorrelation analysis. ICrouzet et al. (I2012D de rived sig- 



nificantly larger errors than Tinetti et al ( 20101). and we 
theref ore find less disagreement with the ICrouzet et aL 
(2012ft spectrum. The discrepancy between ITinetti et al. 
(|2010l) and our result is either due to the intractability of 
NICMOS instrument signatures, or to variability in the 
exoplanetary atmosphere (i .e., clouds at the ter minator). 

If our difference with the ITinetti et all (|2010f) results is 
due to variable clouds at the planet's terminator, it is infor- 
mative to convert the required change in abs orption to an 
equiva lent number of opaque scale heights. ITinetti et al.l 
(|2010l) derive an absorption depth in the spectrum of ap- 
proximately 1150 ppm (from 1.28 to 1.38 /im), whereas we 
measure only ~200 ppm. In terms of equivalent planetary 
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radii, the difference of 950 ppm implies: 

26R P /R P = 0.00095/ '{R 2 p /R 2 s ) (2) 
and adopting R 2 p /R 2 s = 0.017 for XO-lb, we find: 

SR p = 0.0279i? p = 2360 km (3) 

So a ring of height 2360 km would have to be opaque 
with clouds during our measurement, but sufficiently clear 
to allow water absorption at the time of the NICMOS ob- 
servation. The scale height, kT/nmq is about 260 km, 
if we adopt T = 1200X from iMachalek et all (|2008f ) and 
use a molecular hydrogen composi tion (// = 2.3 2 ). In or- 
der to attribute our difference with lTinetti et ail (|2010l) to 
variability of the planetary absorption, requires variable 
clouds around the entire terminator of the planet extend- 
ing over 9 scale heights. We regard this as highly unlikely 
from a meteorological point of view. So we conclude that 
NICMOS exoplanet spectroscopy, as currently analyzed, 
is unreliable. 

As regards HD 209458b, our water ab s orptio n appears 
inconsistent with the results of iBarmanl (|2007f) . Barman 
mentions that his baseline model, that accounts for water 
absorption he identified near 1 /zm, predicts a peak in R p 
of 1.343i?j at 1.4 /im. Estimating the continuum level at 
~ 1.315 Rj from Barman's Figure 1, we project that his 
baseline model would predict a 1.4 /im absorption of about 
580 ppm - about t hree time s what we measure. We re- visit 
the comparison to iBarma n| ((20071) i n Sec. 6.5. 

Using Spitzer transit photometry, iBeaulieu et al.l (|2008| ) 
found evidence for water absorption in HD 209458b. While 
our observations do not overlap the Spitzer bandpasses, 
the IBeaulieu et al.l (|2008l) band-to-band transit depth dif- 
ferences require stronger water absorption than we ob- 
serve, by several times. 

6. INTERPRETATION USING MODEL ATMOSPHERES 

We now address the degree to which our results can 
be affected by stellar activity (star spots, Sec. 6.1), 
then we turn to the interpretation of our results using 
model planetary atmospheres. We begin model interpre- 
tation by implementing a new fast-calculation transmis- 
sion model, and validate it (Sec. 6.2). We also re-analyze 
the HST/STIS optical trans mittance data for HD 209458b 
from lKnutson et"aL1 (|2007al ) (Sec. 6. 3), in order to combine 
those data with our WFC3 results. We then implement 
our new transmittance code (Sec. 6. 4), and then we com- 
pare both the HST/STIS and HST/WFC3 transmittance 
spectra to calculations from other models (Sec. 6.5). 

6.1. Effect of Star Spots 

Prior to discussing our results in terms of exoplanetary 
atmospheric transmission, it is necessary to demonstrate 
that our derived exoplanetary spectra are unlikely to be 
contaminated by stellar activity. Although no distinct 
crossings of star spots are apparent in Figure 3, two ef- 
fects are possible in principle. First, crossings of multiple 
small spots might occur, and could produce a significant 
cumulative effect on the transit spectra. Second, the signa- 
tures of uncrossed star spots - not occulted by the planet 

18 http:/ /solarscience. msfc.nasa.gov/greenwch.shtml 



during transit - could be 'amplified' by the transit phe- 
nomenon. We first consider the latter effect, i.e. possible 
amplification of uncrossed star spots. 

We know from solar observations that the cool um- 
brae of large sunsp ots exhibit water vapor absorption 
( Wa llace et al.l fl995l) . So the spectra of solar- type stars 
will exhibit 1.4 /im water absorption, albeit at a very low 
level. When a planet transits without crossing spot um- 
brae, it increases the relative fraction of the unocculted 
disk that contains spot umbrae, hence it increases the rel- 
ative depth of stellar water absorption during transit. It 
can be shown that the magnitude of the effect masquerad- 
ing as exoplanetary transit spectra is closely approximated 
as A5e, where A is the fractional coverage of the stellar disk 
by umbrae, S is the relative absorption depth of the water 
band in the umbral spectrum, and e is the transit depth. 

We first estimate the magnitude of 5 in the above expres- 
sion. Ideally, we would measure S from observed spectra 
of sunspots at 1.4 ^m, but we do not possess such spec- 
tra due to poor telluric tran smittance. From models of 
cool stars (|Allard et al.ll200"ol ). we estimate that sunspot 
umbrae would exhibit about 30% relative absorption (line 
core to continuum) in the 1.4 /im water band. We check 
that estimate from other properties of sunspots. From 
umbral sp ectra measured in the red -optical continuum 
(U87^mV IPenn fc MacDonaldl (|2007ft found that umbrae 
are about 0.39 as intense as the surrounding photosphere, 
averaged over the solar cycle (umbrae are slightly darker 
at solar maximum, about 0.35 of the photosphere). Con- 
verting this ratio to a brightness temperature, we obtain 
T c = 4380K for the umbral continuum. We estimate the 
temperature in the water line-forming region as T water — 
3200K, from the molecular rotational temperature of th e 
water lines (|Wallace et al.lll995l : iTereszchuk et aLI l2002f ) . 
Applying those brightness temperatures to 1.4 /im, we es- 
timate that the umbral IR continuum intensity is about 
0.52 of the photosphere, and the water band core would 
be about 0.21 of the photospheric intensity if it were opti- 
cally thick. That is about a factor of t wo stronger than w e 
estimated from the cool star models (|Allard et al.l [2000). 
which is reasonable. We conservatively adopt the greater 
(i.e., optically thick) value. The amplitude of the water 
band in the umbral spectrum is about 6 — 0.31 (0.52-0.21) 
of the photospheric intensity. 

HD 209458 and XO-1 are indica ted to have ave r age ac - 
tivity levels in the compilation of iKnutson et al.l (|2010f ). 
About half of solar-type stars in the Kepler sample are 
more active than the S un based on variations in optical 
light ((1 asri et ail 120101 ). We therefore conclude that the 
Sun is a reasonable analogue for HD 209458b and XO-1, 
and we calculate the average value for coverage of the so- 
lar disk by spot umbrae. We use the monthly compila- 
tion of sunspot areas since 1874 given by Marshall Space 
Flight CenteiP^I. Averaging these data, we find that the 
fractional solar disk coverage by sunspots (including their 
penumbrae) is 268 ppm on average. At the maximum of 
the solar activity cycle the typical coverage value is about 
3000 ppm, more than an order of magnitude larger. About 
32% of t he spot area is due to umbrae, based on sunspot 
studies (iBrandt et al.l [l990l) . Therefore the umbral disk 
fraction A is typically 82 ppm, increasing to 960 ppm at 
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the strongest solar maxima. Umbral areas can be diffi- 
cult to estimate, due to scattered l ight and limited spa- 
tial resolution. The seminal work of iHoward et al.l (pL9S4) 
found a significantly lower (300 ppm) disk coverage at so- 
lar maximum based on Mt. Wilson photographic plates 
(1921-1982), but about the same average coverage that we 
infer here. 

For HD 209458b. we have e = 0.0146 (Knutson ct al. 
[2007ah . so ASe is 0.4 to 4.3 ppm, and the numbers for 
XO-1 are similarly small. Even if these stars are several 
times more active than the Sun at solar maximum, the 
amplification-of-star-spots effect is not a major contribu- 
tor to our measured water band depths. 

We now turn to the possible cumulative effect of crossing 
small star spots during transit. Occultation of water ab- 
sorption in star spots would produce an apparent emission 
in the derived exoplanetary spectrum, or it would weaken 
real exoplanetary absorption. The worst case effect would 
occur for a star at the peak of its activity cycle, when all 
of the spots on the stellar disk were occulted by the planet 
during the specific phases of our WFC3 data (Figure 3). 
The star spots would have to be distributed so that they 
were all crossed during our partial coverage of the transit, 
and their sizes would have to be sufficiently equal to pre- 
vent a single large spot-crossing being visible. In this un- 
likely case, the planetary absorption could be weakened by 
0.31(960)= 298 ppm, a significant effect. However, if 960 
ppm of star spot umbrae were occulted during transit, then 
the white-light transit in Figure 3 would be more shallow 
because the umbral continuum is fainter than the photo- 
sphere. That effect would be approximately 960(0.52)= 
500 ppm in transit depth. It would decrease our value of 
R P IR S by approximately 0.002 for both HD 209458b and 
XO-lb. Those decreases are unlikely, give n our agreement 
with Spitzer IR transits of HD 209458b (|Crossfield et al.l 
2D 1 2 ■■ and with an in dependent transit of XO-lb analyzed 
bv lBurke et~aTl (|2010l ) (Sec. 4.1). Moreover, if the weak ab- 
sorptions we derive for these exoplanets are due to mask- 
ing by star spots of intrinsically stronger exoplanetary ab- 
sorption, then the specific and unlikely spot-occultation 
circumstances described above would have to apply inde- 
pendently to both stars. 

We conclude that neither occulted nor unocculted star 
spots have significantly affected our results, and we are 
measuring real exoplanetary water in transmission. 

6.2. Validation of a Spectral Transmittance Code 

In Sec. 6.3, we interpret our results using a new spec- 
tral transmittance code. This code is intended for rapid 
line-by-line calculation of transmittance spectra, so as to 
explore parameter space when varying mixing ratios, cloud 
heights, atmospheric temperature, etc. We intend to use 
it for future WFC3 investigations as well as to illuminate 
the present results. Here, we describe the code, and the 
tests we have conducted to validate it. 

Our transmitt a nce code is based on the work of 
Ric hardson et al.l (|2003T ); consequently many of the algo- 
rithms it uses were tested previously. We have modified 
the original code to use the slant-path geometry appro- 
priate for transmittance spectra. The code uses a layered 
hydrostatic equilibrium atmosphere, with pressures spaced 
equally in the log from 1 to 10 7 dynes cm -2 (adjustable). 



It incorporates co ntinuous opacity due to collision-induce d 
absorption by H 2 (jZheng fe Borvsowll995tlBorvsowl2002l ). 
and it includes a provision for calculating number densities 
of major molecules in thermal equilibrium. Nevertheless, 
we here specify a depth-independent mixing ratio of water 
in an ad-hoc fashion for the purpose of exploring param- 
eter space. We i nclude opacity for water usin g the line- 
by-line data from lPartridge fc Schwenkel (|1997|) . To speed 
the calculations, we sum and pre-tabulate the strengths 
of the extremely numerous water lines within wavelength 
bins, and represent each bin by a single line having the 
strength of the total. We choose a bin width of 0.1 cm -1 
(2 x 10~ 5 ^m at 1.4 fim), much smaller than our WFC3 
resolution. The average within each bin is represented by 
a Voigt profile with a damping coefficient of 0.1 cm -1 per 
atmosphere. This averaging is valid only at a single tem- 
perature, because each bin in wavelength will contain lines 
of different lower state excitation, which adds a non-linear 
factor. We thus pre-tabulate a separate line list for each 
isothermal model. From the line and continuous opacity, 
we numerically integrate along the un-refracted path that 
passes tangent to each layer, to calculate the optical depth 
r and extinction exp(-r) over that path. We adopt a 
source function of zero, i.e., we neglect self-emission by the 
planetary atmosphere. From the optical depths, we calcu- 
late the effective blocking area of the planet as a function 
of wavelength. 

To validate this code, we conducted two tests. First, 
we specified an ad-hoc continuous opacity proportional to 
A~ 4 , a nd we conducted a test described by IShabram et al.l 
(|201lh . This test involves verifying the analytic relation 
dRp/dln A = — 4H, where R p is the calculated wavelength- 
dependent transit radius, and H is the pressure scale 
height. Using a 1000-layer isothermal atmosphere with a 
constant scale height, we calculated the slope of Rp versus 
In A, and we find agreement with the slope of -4 to within 
0.58%. For atmospheres having 500 and 200 layers, our 
precision on the slope is 1.28 and 3.42%, respectively. Our 
grid of spectral calculations to interpret our WFC3 spectra 
(Sec. 6.3) uses the 200-layer version, which is more than 
adequate for the present purpose, as our second test now 
establishes. 

For the second test, we compared our calculated wa- 
ter transmission spectrum to an independent calculation 
from J. J. Fortney. The Fortney mode l is ne w, but uses 
the methods described in lFortnev et al.l (|2010t ). This com- 
parison used our 200-layer version of the code, and an 
isothermal model having T = 1500K, a surface gravity of 
10 meters sec -2 , and a water mixing ratio (independent of 
pressure) of 0.00045. The Fortney calculation uses differ- 
ent algorithms, different layering, and different numerical 
approximations to represent the same atmosphere. Results 
from the two calculations are in excellent agreement, and 
shown in Figure 10, where the (smoothed) results from the 
new code track the Fortney calculation very closely. As an 
additional test, we repeated the Figure 10 comparison us- 
ing a calculation by N. Madhusudhan, adopting different 
atmospheric parameters, and we again achieved excellent 
agreement (not illustrated here). We therefore conclude 
that our new transmittance code is validated for the pur- 
pose of this paper. 
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6.3. Re-analysis of the HST/STIS Photometry 

We want to utili ze our results together w ith transmission 
spect ra derived by IKnutson et all d2007al) from HST/STIS 
data. iKnutson et al.l ( 2007al ) give errors for R p and R s sep- 
arately, but a large fraction of those errors are common- 
mode between R p and R s . To clarify the error on R p /R s 
from the STIS data, we were motivated to re-analyz e that 
photometry (from Table 1 of IKnutson et al. 2007al). Our 
re-analysis adopts many of the Knutson et al.l ~ 2007al ) pa- 
rameters without attempting to vary them. Specifically, 
we adopt their transit center times, their non-linear limb 
darkening coefficients, and their band-to-band differences 
in R p /R s . We fit their photometry for all 10 bands simul- 
taneously, using a Markov Ch ain Monte C arlo (MCMC) 
method with Gibbs sampling (|Fordl [20051 ) . to determine 
the band-averaged value of R p /R s , as well as a/R s and 
i. We in itialize our MCMC ch ains using the parame- 
ters from IKnutson et al.l (|2007al ). and these good start- 
ing values allow our chains to converge rapidly (in ap- 
proximately 1200 steps). However, the chains run slowly, 
because we are fitting to 10 bands simultaneously using 
nonlinear limb darkening. We therefore run 4 indepen- 
dent chains simultaneously, each having 60,000 steps, and 
we combine their posterior distributions. The combined 
distributions are closely approximated by Gaussians, and 
yield R p /R s = 0.1210 ± 0.0001. The small error in R p /R s 
reflects the large volume of data over 10 wavelength bands. 

Our re-determined value for the band-aver a ged R p /R s 
is in mild disagreement with IKnutson et al.l (|2007a), in 
the sense that our result favors a larger R p by 1.5 times 
their error for R p . However, our re-determined STIS value 
for R p /R s is in good agreement with our WFC3 value, so 
our analysis is self-consistent. As noted in Sec . 4.1, we also 
agree with other IR determin ations of R p / R s (jBurke et alj 
2010; Crossfi eld et al.l [20121 ), except for the values from 
Beaulieu et alj (|2008l) . 



6.4. Interpretation of HD 209458b Using Model 
Atmospheres 

We now turn to exploring what these new observations 
imply about the atmospheres of these planets. In the case 
of XO-lb, the scaling factor required to fit the Spitzer- 
based model to the data on Figure 8 was 0.84, versus 0.57 
for HD 209458b. These factors do not per se indicate dis- 
crepa ncies with the conclu s ions from the Spitzer investiga- 
tions (jKnutson et al.ll2008t iMachalek et al.ll2008ft . because 
molecular absorption during transit is much more sensi- 
tive to conditions such as clouds and haze. Given that 
the scale factor for XO-lb is close to unity, and consid- 
ering the larger noise level of that spectrum compared to 
HD 209458b, we conclude that XO-lb is a dequately de- 
scribe d by extant models, not counting the iTinetti et al.l 
(2010) calculations. In contrast, the relatively small scale 
factor required for HD 209458b, combined with the lower 
noise level of those data, motivate us to inquire further 
what our results imply for the atmosphere of that planet. 

In principle, it would be possible to analyze 
our water transmission spectrum s imultaneously with 
sodium absorption measurements (|Charbonneau et al.l 
l200llSing et al.ll2008t ISnellen et al.ll2009h . CO absorption 
data ( Snellen et al.ll2010D an d Spit zer secondary eclipse 
photometry ( Knutson et alj 120081 ) and spectroscopy 



ijRichardson et all 120071 : ISwain et al.l l2008bl ). Such an 
analysis could incorporate guidance from hydrodynamic 
models dShowman et all l2008t iDobbs-Dixon et all 120121: 
IShowman et alj 120131 : iRauscher fc Menoul l2013f ) to ac- 
count for the longitudinal transfer of stellar irradiance, 
and could explore the full p arameter space of composition 
and t emperature structure (|Madhusudhan fc Seagerl l2009. 
2010j), as well as the effect of clo uds and hazes at the ter- 
minator (jParmentier et afll2013f ). However, we do not at- 
tempt such an ambitious investigation here. Instead, we 
compare the combination of HST/STIS and HST/WFC3 
measurements to theoretical transit spectra from various 
models without much fine-tuning, and we discuss the dis- 
crepancies so as to guide the general direction of more 
exhaustive atmospheric modeling in the future. 

The weakness of the water band we observe is reminis- 
cent of the first detec tion of this (indeed, o f any) exoplan- 
etary atmosphere by ICharbonneau et"aU (|2002h who ob- 
served the sodium 'D' line doublet near 0.5893 ^im using 
HST/STIS. To account for t he we akness of the sodium 
feature, ICharbonneau et "all (j2002) varied the height of 
an opaque c l oud layer, using the model described by 
iBrown et alj (|2001f h We here perform a similar calcula- 
tion, using our spectral transmittance code described in 
Sec. 6.2. We vary both the pressure level for the top 
of an opaque cloud layer, and the mixing ratio of water 
vapor, and we measure the strength of the 1.4 /im band 
at about 300 points over a 2-D grid, calculating a full 
transmittance spectrum (as per Figure 10) at each grid 
point. The question of the temperature structure at the 
terminator of the planet - needed for this calculation - 
is problematic, i n spite of significant work on this topic 
( Sin g et al~ll2008[ ). Our results are not very sensitive to the 
T(P) relation, so we prefer to use a simple isothermal at- 
mosphere, ba sed on the observed b rightness temperatures 
from Spitzer (Crossfiel d et al.l l 20 12). Ideally, we would av- 
erage the day-side and night-side brightness temperatures 
from Spitzer to arrive at an estimate of the temperature at 
the terminator. Although Spitzer around-the-orbit obser- 
vations for HD 209458b have been obtained, they are still 
under analysis. Therefore we adopt the same relative day - 
to-night change as for HD 189733b (jKnutson et al.ll2007bl) . 
and apply that relati ve variation to the HD 209458b day 
side temperature from [Crossfield et al. ( 2012). This yields 
an estimate of 1200K for the terminator, and Figure 11 
shows a grid of band amplitudes calculated at that tem- 
perature. 

The blue line with ±lcr error limits on Figure 11 is the 
amplitude of the observed band. Since the models have 
fine-scale structure with wavelength, and the observations 
have point-to-point noise, we define a band amplitude for 
Figure 11 by averaging over the same wavelength regions 
mentioned in Sec. 4.4 and used for Figure 6, i.e., 1.27-1.30 
and 1.36-1. 44 /im. This results in an observed band am- 
plitude of 176±25 ppm. Note that the amplitude of the 
absorption near the actual bandhead at 1.38 /jm is slightly 
higher, about 200 ppm. 

Contours of constant mixing ratio are included on Fig- 
ure 11. The solar- abundance contour (drawn in red) inter- 
sects the observed band amplitude only where the cloud 
top pressure is low, near 1.5 mbars. T hat is similar to, but 
less extreme, than the sodium case; ICharbonneau et al.l 
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(2002) remarked that if the sodium weakness is attributed 
solely to clo uds, then it 'w o uld re quire.. .clouds tops above 
0.4 mbar.' iFortnev et al.l (|2003f) concluded that silicate 
and iron clouds could reside at pressures of several mil- 
libars in the atmosphere of HD 209458b. As concerns al- 
ternate explanations, note that Figure 11 implies a semi- 
forbidden region, where no contours pass into the lower 
right of the plot. The cloud-top pressure levels depicted 
on the right of Figure 11 imply a clear atmosphere, and 
even our lowest modeled mixing ratio (-5.2 in the log) is 
not sufficient to weaken the band to account for our ob- 
servations if the atmosphere is clear. The total column 
density along the line of sight at high pressures in the tan- 
gent geometry is so large that even unrealistically small 
mixing ratios are insufficient to weaken the band to the 
observed degree, in the absence of other water-destruction 
mechanisms such as photolysis (unlikely in the deep at- 
mosphere). We conclude that we are not observing a clear 
atmosphere. 

Figure 12 shows our full HD 209458b transit depth spec- 
trum (i? 2 /i? 2 ), combining both our WFC3 and re-analyzed 
STIS results. The combination of these data span wave- 
lengths from 0.2 to 1.6 /xm with a consistent observed lower 
envelope and overall level in i? 2 /i? 2 . We first add the 
caveat that systematic differences might still remain be- 
tween the overall level of the STIS and WFC3 transit 
depths, in spite of the seeming consistency. Nevertheless, 
Figure 12 represents the best composite optical/near-IR 
transmission spectrum of HD 209458b to date, so we pro- 
ceed to ask what it reveals about the exoplanetary atmo- 
sphere. 

Now, we compare this combined spectrum to two mod- 
els. First, we used a grid of spectra by Adam Bur- 
rows, based on a 1200K isothermal temperature structure, 
as used a bove. The grid utili zes t he methodolog i es de- 
scri bed by Burrows et all (1200 lh and lBurrows et all ((2010) 
and iHowe & Burrows! (|2012j Tbut it incorporates differem, 
amounts of extra gray opacity. We interpolate in this grid 
to find that an extra opacity of 0.012 cm 2 g _1 matches 
the 1.4 /xm water absorption at the bandhead, and pro- 
vides suitably low absorption at 1.15 /im. The lowering of 
the 1.15 /im absorption occurs because that intrinsically 
weaker band requires a longer path length to produce sig- 
nificant absorption, and long path lengths are masked by 
the extra opacity. A Burrows isothermal model having no 
extra opacity (not illustrated) shows a much more promi- 
nent peak at 1.15 /im. 

The profile of the 1.4 /jm band is not matched optimally 
by the isothermal models, not as well as on Figure 8 for 
example. The real absorption line of sight passes through 
different temperatures on day and night hemispheres of the 
planet. Figure 8 accounts (crudely) for different tempera- 
tures along the line of sight, not included in the isothermal 
model for Figure 12. Including that line-of-sight tempera- 
ture variation may be essential to matching the band pro- 
file. 

The Burrows model on Figure 12 is sufficiently 
high-resolution in wavelength to permit meaning- 
ful comp arison with the s odium abs orption mea- 
sured by ICharbonneau et ahl (|2002h and ISnellen et al.l 
(I2009D. We plot the ' narrow' band absorption from 
ICharbonneau et"al] ((2002) on Figure 12 (triangle point 



with error bar). The lSnellen et al.l (I2009D resu lts (not plot- 
ted) are consistent with ICharbonneau et al.l (|2002h . con- 
sidering the different bandpasses. Integrating the Burrows 
model h aving 0.012 cm 2 g -1 extra opacity over the band 
used by ICharbonneau et all ((2002) (red diamond point) 
produces agreement within the error bar. 

One aspect of the observations that are not reproduced 
by the simple isothermal Burrows model with gray opac- 
ity, is the tendency toward increasing radius in the blue 
and UV, at the left edge of Figure 12. An increase 
of transit radius at short wavelength may be related to 
the abs orber that causes a temperature inversion in this 
planet (jBurrows et al.ll2007t ). It may also be produced by 
Rayleigh scattering from a population of small particles, 
that we do not include in the Burrows calculations for Fig- 
ure 12 (but Rayleigh scattering by molecules is included). 
In this regard, we overplot a model from Ian Dobbs-Dixon 
(b lue line on Figure 1 2 ), bas ed on the methods described 
in iDobbs-Dixon et ail ( 2012f) . This model uses a full ra- 
diative hydrodynamic treatment of the temperature struc- 
ture, which may explain why it produces a better (but 
not perfect) account of the 1.4/xm band profile. It has 
no extra gray opacity, but it incorporates extra opacity 
of 0.004cm 2 g _1 at 0.8 /im, with a A -4 dependence. Be- 
cause that Rayleigh opacity is concentrated at short wave- 
lengths, we must scale the modulation in the modeled 
spectrum downward by a factor of 3 to match the ob- 
served 1.4 /im band. That scaling is unphysical, but it 
allows us to judge the relative importance of gray versus 
Rayleigh opacity that will be needed to match the obser- 
vations. After scaling, the blue line produces relatively 
good agreement with the 1.4 /im data, but overestimates 
the 1.15 /im feature as well as slightly overestimating the 
increase in absorption in the blue and UV. 

Finally, we point out one notable discrepancy in the 
model comparisons. The STIS point near 0.95 /zm can- 
not be reproduced by models while s till being consis tent 
with our 1.4 /im band measurement. iBarma 3 (pool ar- 
gued for water vapor in HD 209458 b based in part on the 
STIS data near 0.95 /mi, but iKnutson et al.l (|2007aD did 
not claim water detection from those same data. 

From the above comparisons, we conclude that: 

• A uniformly distributed extra opacity, gray in 
wavelength dependence, of approximate magnitude 
0.012 cm 2 g _1 is needed in isothermal models at 
1200K in order to match our observed water transit 
absorption for HD 209458b at 1.4 /im, and to pro- 
duce the weakness of water at 1.15 /zm, and account 
for weak sodium absorption in the optical. 

• Models that include realistic temperature distribu- 
tions along the line of sight may be required to 
match the band profiles of the water absorption. 

• Our results are consiste n t wit h the situation de- 
scribed by iPont et al.l (|2012l ) for HD 189733b, 
wherein weak molecular absorptions are superposed 
on a spectrum whose broad variations with wave- 
length in transit are dominated by haze and/or 
dust opacity. However, the extra opacity required 
for HD 209458b needs to be grayer than the strong 
Rayleigh component needed for HD 189733b. 
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7. SUMMARY AND FURTHER IMPLICATIONS 

We have demonstrated that WFC3 spectroscopy using 
the new spatial scan mode can yield exoplanetary spectra 
whose error level falls significantly below the 0.01% level. 
We detect 1.4 /zm water in both HD 209458b and XO-lb 
at a relatively low level of absorption, only 200 ppm at 
the 1.38 /zm bandhead. Our results for XO-lb contradict 
the much larger absorption derived by Tinetti et al.(2010), 
and we conclude that NICMOS spectr oscopy as currently 
analyzed is unreliable. iFortnevI (|2005l ) predicted that gi- 
ant exoplanetary atmospheres would contain condensates 
and hazes that would 'lead to weaker than expected or un- 
detected absorption features'. For HD 209458b, an atmo- 
sphere with 0.012 cm 2 g -1 of extra opacity, gray in charac- 
ter, is required in order to match the subtle water absorp- 
tion we detect an d the sodium absorption in the optical. 
iPont et ail (|2012[ ) have argued that molecular absorptions 
in HD 189733b are weakened by haze and/or dust opacity, 
and our results suggest a similar situation for HD 209458b. 
One difference is that haze and/or dust opacity is grayer 



for HD 209458b, with a weaker Rayleigh component as 
compared to HD 189733b. 

The capability to derive transmittance spectra for plan- 
ets transiting bright stars, at an error level below 0.01% 
in a single HST visit, opens new scientific horizons. For 
example, it becomes feasible to monitor meteorological 
variability of conditions at the terminator of the planet. 
The transmittance spectra of su per-Earths should be de- 
tectable using multipl e visits (|Bean et"all l2010l 120111 : 
iKriedeberg et all I2013D . even for high molecular weight 
atmospheres. Applying the spatial scan to secondary 
eclipses, it should be possible to confirm the existence of 
atmospheric temperature inversions via low-noise detec- 
tion of the water band profile in emission rather than in 
absorption. Finally, we note that there is much discussion 
in the community concerning dedicated space missions to 
characterize exoplanetary atmospheres using transits. The 
design of such missions should prudently consider that 
molecular absorptions may be considerably weaker than 
are modeled using clear atmospheres. 
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APPENDIX 



PRODUCTION OF 2-D SPECTRAL FRAMES 

We here describe how the spectral frames shown in Figure 1 are constructed from the sample-up-the ramp data frames, 
and how discrepant pixels and energetic particle hits are corrected. Each non-destructive sample of the detector is provided 
to observers as an extension in a FITS file (the *ima.fits files available from the Multi-Mission Archive at Space Telescope, 
MAST). Normally the *ima.fits files are processed by an analysis pipeline at the Space Telescope Science Institute (STScI), 
which fits the slope and makes the resulting intensities available in *flt.fits files. However, in the case where the source 
is rapidly moving, as it is for the spatial scan mode, the normal analysis pipeline is inapplicable because a given pixel 
is not always viewing the same celestial scene, and that invalidates the pipeline's fit to the slope at each pixel. In scan 
mode, there are two alternative methods that can be used for initial data processing. One method is to subtract the 
first read of the detector from the last read, thus producing an image of the total accumulated electrons on the detector 
during the spatial scan. However, we use an alternate method that offers several advantages over a simple 'last minus 
first' subtraction. 

For a given exposure, let M be the number of times that each pixel is sampled, and let i = 1, M index the individual 
samples. For each i < M — 1, we subtract sample i from sample i + 1, thus forming differences Di = Cj+i — C'i where C 
denotes the charge on a given pixel in electrons. Because scanning occurs parallel to columns, these differences reveal the 
spectrum of the star over a limited range of rows on the detector. Considering the difference Di, we zero the rows not 
containing the target star, using a top-hat mask whose width extends to the wings of the stellar point spread function 
(PSF) - typically 15 pixels (1 pixel =0.14 arcsec). We then form the sum: 



where j = 1, N indexes the N times that the scan is repeated over the duration of the transit. Each sum is the image 
of the spatial scan that we use in further analysis. 

The above methodology reduces the effect of sky background. Because of the masking applied to each of the Di, sky 
background not immediately adjacent to or underlying the star is zero-ed, and does not accumulate in the summation 
in Eq. (1). Only the background near the star in each Di survives this procedure. However, the background flux that 
accumulates in the brief interval between consecutive reads of the detector is small (< 0.1%) compared to the fluxes from 
these bright stars, and is not removed at this stage of the analysis. Instead, it is removed by our spectral derivation 
procedure, described in Sec. 4. Advantages of the above procedure are that it minimizes the effect of hot pixels and 
energetic particle events that would otherwise overlap the scan. Moreover, it allows for discrimination against other stars 
that are spatially resolved, but would overlap in a simple last-minus first difference. 

Although the procedure described above minimizes the effects of discrepant pixels, it does not completely eliminate 
them. Normally, hot and transient pixels are identified and corrected via a numerical median filter applied to the time 
history of each pixel. The spatial scan rate variations complicate that procedure because they contribute to the intensity 
fluctuations of every pixel versus time, and could interfere with the median filter process. We therefore apply a 5-point 
median filter to the ratio of a given pixel intensity to the total intensity in that row of the detector. This ratio cancels 
the scan rate variations (Figure 1), and isolates the behavior of the pixel itself. We apply the median filter in a two-pass 
process. The first pass corrects pixels that are discrepant by more than 10cr from the median value, where a is the 
standard deviation of the difference between the time history of the pixel and the median-filtered version of that time 
history. The first 10a pass serves to eliminate very large fluctuations that might perturb the calculation of a. The second 
pass uses a lower threshold, correcting pixels that are discrepant by more than 3a. For HD 209458b, 0.15% of the pixels 
are corrected by this procedure, and 0.04% in the case of XO-1. 



M-l 




(1) 



i=l 



-0.01 0.00 0.01 

Fractional difference 

Fig. 1. — Images of the spatial scans at representative times. Left to right, the images are: a scan of XO-lb (20 pixels scan height), a scan of HD 209458b (228 pixels scan height), and 
the difference between two consecutive scans of HD 209458b. The difference images shows a striped appearance due to small variations in the scan rate under control of the Hubble fine 
guidance system. The plot shows the fractional intensity fluctuations as a function of scan position (Y-axis) for the difference image. 
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Fig. 2. — Example of the most prominent instrument-related systematic effect in staring-mode data from our Cycle-18 WFC3 program. 
This example zooms-in on the egress portion of a transit of WASP- 18b, using the integral intensity over each grism spectrum. The space 
between groups of observed points is due to the time needed to transfer the data buffer. Within each group of points, the intensity increases 
in a pattern shaped somewhat like T'. We call this pattern the 'hook'. 
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Fig. 3. — Transits observed in spatial scan mode: HD 209458b (top) and XO-lb (bottom), integrating over the entire grism bandpass ('white 
light'). No attempt was made to remove systematic effects by divide-oot methodology I Berta et al these are purely 'as observed'. The 

blue cur ves are fit to the data by v arying R v /R s and a correction to the time of central transit, but fixing other parameters at the values 
given by iKnutson et al] l|2007al ) and iBurke et alj 1(20101) . 



18 



Deming et al. 



Column 

20 40 60 80 100 120 



6x10" 



on 

a 
o 

§ 4xl0 5 

<u 



2x10" 



1.002 

o 

1.000 

Pi 

0.998 



0.996 



T~ 





1.1 1.2 1.3 1.4 1.5 
Wavelength (microns) 



1.6 



1.7 



Fig. 4. — Top panel: example of the grism spectrum of XO-1. The black line is a spectrum at a single time, showing the roll-off in grism 
response shortward of 1.1 fim and longward of 1.65 (im. The red curve is the out-of-transit template spectrum (see text), shifted upward by 
1% for clarity of illustration. These spectra have been smoothed using a Gaussian kernel having FWHM = 4 pixels. Lower panel: difference 
between the single smoothed spectrum and best-fit template, normalized by the intensity in the template spectrum. The vertical dashed lines 
define the wavelength range used in the analysis of XO-lb. (A similar range was used for HD 209458b.) 
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Fig. 5. — Example of undersampling in the WFC3 spectra. Shown are zoomed-in portions of two spectra of HD 209458, separated by 3.1 
hours of time (in orbits 1 and 3 of visit 26 in our program). The dip at 1.28 /an is the Paschcn-beta line in the star. In the earlier spectrum 
(solid line), the line core appears flattened because the line is positioned mid- way between two columns of the detector. In the later spectrum 
(dashed line), the line core is sharp because the different overall shift in wavelength places it centered on a column. 
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Fig. 6. — Top panel: differential transit at the peak wavelength of the exoplanetary water absorption in HD 209458b, by averaging residuals 
from the template spectrum fit (as per Figure 4) over the wavelength range 1.36-1.44 fim. Note the residual deeper-than-avcrage transit depth. 
Bottom panel: the same as the top panel, but for the wavelength range 1.27-1.30 /^m where little water absorption is expected. In this case 
the differential transit depth is negative, i.e. less deep than the average over the total wavelength range analyzed. 
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Fig. 7. — Error analysis for our derived exoplanetary transmission spectra. Each panel plots the standard deviation of the observed noise 
in our differential transits, after removing the best-fit amplitude. The noise is shown as a function of bin size. The blue lines are the relations 
expected for photon noise based on the number of detected electrons, and accounting for the effect of smoothing the grism spectra (Sec. 4.3). 
The blue lines have a slope of -0.5 due to the expected inverse square-root dependence of the noise; the measured points are in good agreement 
with that expectation. 
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Fig. 8. — Our results for transmission spectra for HD 209458b and XO-lb in the WFC3 bandpass, compared to model predictions (blue 
lines). The spectral resolving power of these measurements is A/<5(A) ~ 70. The amplitude of the 1.4 (im water absorption is about 200 
parts-per-million (ppm) in both cases, but the errors are smaller for HD 209458b due to the greater photon flux. 
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Fig . 9. — Comparison o f our transmission sp ectrum of XO-lb (solid points, see Figure 8) with the NICMOS results from iTinetti et al, 
12010) (black squares) and ICrouzet et al] 1 2012 ') (blue squares). (Our data have been offset in transit depth for clarity.) Our WFC3 water 
absorption is of much smaller amplitude than seen in the NICMOS data (see text for discussion). 
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Fig. 10. — Validation of our transmission spectral model (blue line) versus an independent calculation from lFortnev et alj 12010) for an 
isothermal model at 1500K. 
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Fig. 11. — Depth of the 1.4 /im water band versus cloud top pressure for an isothermal model at 1200K. Lines represent different mixing 
ratios of water. From top to bottom the log of the mixing ratios vary from -2.0 to -5.2 in increments of -0.2. The red contour is the mixing 
ratio expected for solar abundance. Small black points are the actual calculations from our transmittance model. The blue line is our observed 
band depth of 176 ppm, with ±1ct errors (dashed lines). For this Figure, band depth is defined as the average transit depth from 1.36-1.44 /an 
minus the average depth from 1.27-1.30 fim. 
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Fig. 12. — Transmission spectrum of HD 209458b d erived from Hubble spe ctroscopy. Our WFC3 results are the solid points. The open 



squares are our reanalysis of the STIS bands defined by 



Knutson et al 



Charbonncau et al 



1 2007al ) , and the diamond is the narrow sodium band absorption from 



12002). The red line is the transmittance spectrum from an isothermal Burrows model, having an extra opacity of gray 
character and magnitude 0.012 cm 2 g — 1 . The red diamond integrates the red model over the sodium bandpass. The blue line is a Dobbs-Dixon 
model for HD 209458b, with no gray opacity, but with A -4 (Rayleigh) opacity, normalized to magnitude 0.001 cm 2 g -1 at 0.8 /im. Because 
the blue model has no gray opacity, we scale-down the modulation in this spectrum by a factor of 3 for this comparison (see text). 
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Table 1 

Summary of the spatial scan observations 





XO-lb 


HD 209458b 


Time of first scan B JD (TDB) 


2455834.6666 


2456196.0895 


Planetary orbital phase at first scan 


-0.0471 


-0.0566 


Time of last scan BJD (TDB) 


2455834.9419 


2456196.3997 


Planetary orbital phase at last scan 


0.0228 


0.0314 


Number of scans 


128 


125 


Number of HST orbits 


5 


5 


Scan rate (arcsec per sec) [pixels per sec] 


(0.05) [0.41] 


(0.9) [7.44] 


Detectory subarray size 


128x128 


256x256 


Detector reads per scan 


8 


4 


Duration of scan (sec) 


50.4 


32.9 


Signal level on detector (electrons per pixel) 


3.8 x 10 4 


4.4 x 10 4 



Note: planetary orbital phase is denned to be zero at mid-transit. 



Table 2 



Results for Radius Ratios (R p /R s )and Mid-Transit Times 





XO-lb 


HD 209458b 


Mid-Transit Time BJD (TDB) 


2455834.85186 ±0.00017 


2456196.28934 ±0.00018 


Rp/Rs 


0.1328 ±0.0006 


0.1209 ±0.0004 



